Effect of Ceo giant resonance on the photoabsorption of encaged atoms 



(N 

o 

(N 

Q 
CO 



C/2 



o 

CZ5 

>> 

Or 



>: 
o. 

(N 
t— i ; 
ov 

o: 



Zhifan Chen and Alfred Z Msczane 
Department of Physics and CTSPS, Clark Atlanta University, Atlanta, GA 30314 USA 

(Dated: September 10, 2012) 

The absolute differential oscillator strengths (DOS's) for the photoabsorption of the Ne, Ar, and 
Xe atoms encapsulated in the Ceo have been evaluated using the time-dependent-density-functional- 
theory, which solves the quantum Liouvillian equation with the Lanczos chain method. The calcu- 
lations are performed in the energy regions both inside and outside the Ceo giant resonance. The 
photoabsorption spectra of the atoms encaged in the Ceo demonstrate strong oscillations inside the 
energy range of the C60 giant resonance. This type of oscillation cannot be explained by the confine- 
ment resonance, but is due to the energy transfer from the Ceo valence electrons to the photoelectron 
through the intershell coupling. 

PACS numbers: 33.80.-b, 33.80.Eh 



I. INTRODUCTION 

Recently a method to evaluate the absolute differen- 
tial oscillator strengths (DOS's) for the photoabsorption 
of atoms encapsulated inside a fullerene has been de- 
veloped ■ This method takes three steps to calcu- 
late the photoabsorption spectra of the fullerene and the 
cndohcdral fullerene separately. Firstly, the structures 
of the fullerene and endohedral fullerene are optimized 
0]. Secondly, the ground state eigenvalues and eigen- 
vectors are created by solving the Kohn-Sham equation 
self-consistently using supercell [3] and a plane wave ap- 
proach [f|. Thirdly, the linear response of the system 
to the perturbation by an external electric field is de- 
scribed by the quantum Liouvillian equation Q. The 
photoabsorption spectra are evaluated using the time- 
dependent-density-functional-theory (TDDFT) with the 
Liouville-Lanczos approach 0, H| . Finally, the absolute 
DOS's for the photoabsorption of an atom encapsulated 
inside a fullerene can be calculated by subtracting the 
DOS's of the fullerene from the correspondent DOS's of 
the endohedral fullerene at the same photon energy. 

This method has been successfully used to study 
the photoabsorption spectra of the Xe@Cgo [1] and 
Sc3N@C8o [2] molecules. These examples demonstrate 
the great advantage of the method. The method can 
be used to study the atoms located at the center of the 
fullerene [1] and the off-center positions as well [2]. It 
also allows us to evaluate the spectrum in a broad en- 
ergy region. 

In the photoionization studies of endohedral fullerenes 
the most interesting problem is probably the Xe Ad giant 
resonance of the Xe@C6o molecule [9l-fl4j. The theoret- 
ical calculations that used model potentials [To|, EH and 
the TDDFT with jellium approximation [l2[ ignored the 
effect of the Cgo giant resonance because, fortunately the 
Xe 4d giant resonance is found at the energy region far 



from the Ceo resonance. Amusia and Balenkov [151 ] and 
Madjet et al [l6| have studied the Cgo resonance effects 
on the photoionization cross sections of the Xe 5s and Ar 
3p electrons encaged in the Ceo- However, the photoab- 
sorption spectra of the Ar and Xe atoms encapsulated in 



the Cgo is still unknown. In this paper we have evalu- 
ated the absolute DOS's for the photoabsorption of the 
Ne, Ar, and Xe atoms encapsulated in the Cgo in the 
regions both inside and outside the Cqq giant resonance. 
The results demonstrate a strong oscillation and a signif- 
icant increase of the DOS's in the Ceo resonance region 
for all these atoms. 



II. METHOD 

As stated in the Introduction the method of calculation 
used here takes three steps to evaluate the photoabsorp- 
tion spectra of the fullerene and endohedral fullerene. 
In the first step we utilize the DM0I3 software package 
0] to determine the optimized structure of the fullerene. 
Geometry optimization of the Ceo was performed us- 
ing the generalized gradient approximation (GGA) to 
the density- functional-theory (DFT) with Perdew- 
Burke-Ernzerhof (PBE) exchange-correlation functional 
\vi\ along with all electrons double numerical plus po- 
larization (DNP) basis sets and dispersion correction as 
implemented in the DM0I3 package 0. The optimiza- 
tion of atomic positions proceeded until the change in 
energy was less than 5xl0 -4 eV and the forces were less 
than 0.01 eV/A. The optimized structure was then in- 
troduced into a supercell of 18A. The Kohn-Sham equa- 
tion was solved self-consistcntly to create the ground 
state eigenstates and eigenvalues for the total of 240 
electrons and 120 states using the plane wave approach 
0. An ultrasoft pseudopotential using Rappe-Rabe- 
Kaxiras-Joannopoulos (RRKJ) pscudization algorithms 
fl8j , which replace atomic orbitals in the core region with 
smooth nodclcss pscudo-orbitals, has been employed in 
the calculation. The kinetic energy cutoff of 408 eV for 
the wave function and 2448 eV for the densities and po- 
tentials were employed in a standard ground state DFT 
calculation 0. The linear response of the ground state 
to an external perturbation by an electric field was de- 
scribed by the quantum Liouvillian equation 0,0 (atomic 
units are used throughout, unless stated otherwise): 
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i-£ = [H KS (t),p(t)], (1) 

where Hj(s(t) is time dependent Kohn-Sham (KS) 
Hamiltonian and p(t) is one-electron KS density matrix. 
The KS Hamiltonian can be written as 

H KS {t) = -^V 2 + v ext (r, t) + v Hxc {r, t) (2) 

where v ext (r,t) is the external potential and VHxc( T jt) 
is the time dependent Hartree potential plus exchange- 
correlation potential. 

Linearization of Eq. (1) with respect to the external 
perturbation leads to: 

i% = lH% S s,p'(t)] + [v' Hxc (t), Po] + K X M Po] (3) 

where H^f. is time independent ground-state Hamilto- 
nian, v' ext (t) is the perturbing external potential, v' Hxc (t) 
is a linear variation of the Hartree plus exchange- 
correlation potential induced by p'it). p'(t) = p(t) — po- 
The linearized Liouvillian equation is given by 

i^-=L.p'+[vUt),Po], (4) 

The action of the Liouvillian super-operator L on p' is 
given by 

L-p' = [H% s s ,p'] + [v' Hxc [p']{t),p Q ] (5) 

Fourier analyzing Eq. (5) we have: 

(w — L) ■ p'(w) = [v' ext (w),p ] (6) 

if v' ext (r, w) = — E(iu) • r, the response of the dipole to an 
external electric field E(w) is given by 

d i( w ) = a ij E j( w ) ( 7 ) 
The dynamical polarizability, aij(ui) is defined by 

a t3 {u) = -<r l \hlfA_ > (8) 

Eq. (9) indicates that the dynamical polarizability can 
be expressed as an appropriate off-diagonal matrix ele- 
ment of the resolvent of the non-Hcrmitian Liouvillian 
superoperator between two orthogonal vectors [(| . These 
matrix elements are calculated using the Lanczos algo- 
rithm @, Hi- Finally the absolute DOS's of the C 60 are 
obtained from 

After the Ceo calculation a Xe atom was introduced 
into the center of the Cgo- The DOS's of the Xe@C6o was 
evaluated using the same procedure as described above. 
Finally the DOS's of Xe atom encapsulated in the Ceo 
were obtained by subtracting the DOS's of the Ceo from 
the corresponding DOS's of the Xe@Ceo molecule. 



III. RESULTS 

Fig. 1 presents the spectra of a Xe atom encapsu- 
lated inside the Cgo in the energy range of Xe 4d giant 
resonance. Solid curve and dashed curve represent, re- 
spectively the results of the TDDFT [1] and using our 
C 60 model potential [13] • Both the TDDFT and model 
potential results confirm the three main peaks observed 
in the experiment [l4| if the spectra are reduced by a 
factor of 8 or 10 as discussed in Rcfs. [1,13]. 

The photoabsorption spectrum of the Xe atom encap- 
sulated in the Cqo in the energy region of Ceo giant res- 
onance can be evaluated using the same method as men- 
tioned above. It is well known that there are t wo g iant 
resonances in the Cgo photoionization spectrum [19l - [2l| . 
which correspond to a collective oscillation of the 240 de- 
localized valence electrons relative to the ionic cage of the 
Ceo fullerene. These delocalized electrons are distributed 
over the surface but confined in a thickness of a single car- 
bon atom in the radial direction. The valence electrons 
can move coherently in photoexcitation. The strong peak 
around 6.3 eV [l9| is interpreted as collective excitations 
of the 7r electrons, the giant resonance between 15 -25 eV, 
peaked at 22 eV belongs to the a electrons [2(| [2l|. In 
addition to these there are two broad structures around 
respectivel y, 2 9 eV and 32 eV , which are caused by shape 
resonance [2l| . The whole Cgo spectrum can be found 
inRef. [1]. 

Solid and dashed curves in Fig. 2 are, respectively the 
absolute DOS's for the photoabsorption of the encapsu- 
lated and free Xe atom calculated by TDDFT method. 
In the energy range of 65 -140 eV, the solid curve is same 
as that in Fig. 1. Within the Ceo giant resonance re- 
gion the spectrum of the Xe atom encapsulated in the 
Ceo shows strong oscillation. It is noted that the ab- 
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FIG. 1. (Color online) Photoabsorption spectra of a Xe atom 
encapsulated inside Ceo in the energy range of the Xe 4d giant 
resonance. Solid and dashed curves represent, respectively 
the results of our TDDFT [1] and the Ceo model potential 
calculation [13]. 
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FIG. 2. Absolute DOS's for the photoabsorption of the Xe 
atom. Solid and dashed curves represent, respectively the 
TDDFT results of the encapsulated and free Xe atoms. 



FIG. 3. (Color online) Absolute DOS's for the Photoabsorp- 
tion of a free Neon atom. Solid curve and open circles rep- 
resent, respectively our TDDFT calculation and the experi- 
mental data [231 ]. 



solute DOS's for the atom encapsulated in the Ceo was 
obtained by subtracting the DOS's of the Ceo from the 
corresponding results of the Xe@Cgo molecule. There- 
fore if the DOS's of the Cgo fullerene is larger than the 
DOS's of the Xe@C6o molecule the curve has negative 
values. This means that the Xe atom disturbed the Ceo 
resonance and reduced the DOS's of the Ceo- Fig. 2 
indicates that the negative value usually happens in the 
low energy range. If the photon energy cannot ionize the 
Xe hp electron the Ceo valence electrons hardly transfer 
their resonance energy to the Xe electrons through inter- 
shell coupling. Otherwise the energy transfer from Ceo 
valence electrons to the photoelectron is strong. Because 
of this, the peak at 25 eV may be related to the Xe 5s 
ionization (see Fig. 2). The ionization threshold for the 
5s of a free Xe is 25.4 eV. 

Strong oscillations as in Fig. 2 can also be found in the 
spectra of the Ne and Ar atoms encapsulated in the Ceo- 
In the following, we first evaluate the DOS's for the pho- 
toabsorption of the free Ne and Ar atoms. The ultrasoft 
pseudopotentials for these atoms have been created using 
the PBE [r7[ exchange-correlation functional and includ- 
ing the relativistic effect and nonlinear core corrections 
[221 . The valence electrons for Ne and Ar arc respec- 
tively, 2s, 2p and 3s, 3p. Solid curves in Figs. 3 and 4 
are, respectively the calculated absolute DOS's for the 
free Ne and Ar atoms. The open circles in Figs. 3 and 4 
represent the correspondent experimental data 23j. The 
measurements used the low-resolution dipole spectrome- 
ter. The data started from the first ionization potentials, 
21.6 eV of the Ne 2p state and 16.0 eV of the Ar 3p state. 
The present TDDFT calculations are in good a gree ment 
with the experimental results reported by Ref. [23| . 

The solid curves in both Figs. 3 and 4 rise rapidly 
initially. The experimental spectra of Ne and Ar, reach 
their maxima, respectively at the 31.5 and 21.6 eV. Af- 
ter reaching its maximum the Ar curve drops quickly to 




Photon energy (eV) 



FIG. 4. (Color online) The symbols have the same meaning 
as in figure 3, except that the curve and open circles represent 
the free Ar atom. 



approach the Cooper minimum around 47 eV. The 3s 
ionization threshold of a free Ar atom is about 34.7 eV. 
At this energy a small peak is found in the theoretical 
curve. 

The curves in Figs. 3 and 4 demonstrate that our ultra- 
soft pseudopotentials describe the photoabsorption pro- 
cesses of the Ne and Ar atoms very well. The photoab- 
sorption spectra of the endohedral fullerenes, Ne@Cgo 
and Ar@C6o have also been calculated using the same 
TDDFT method. The DOS's of the Ne and Ar atoms 
encapsulated in the Cgo are obtained by subtracting the 
DOS's of the Ceo from the correspondent DOS's of the 
endohedral fullerenes. 

Dashed and solid curves in Figs. 5 and 6 represent, 
respectively the DOS's for the photoabsorption of the free 
and encaged Ne and Ar atoms. The solid curves in Figs. 
5 and 6 demonstrate the strong oscillations in the region 
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FIG. 5. (Color online) Solid and dashed curves are respec- 
tively, the calculated absolute DOS's for the photoabsorption 
of the encapsulated and free Ne atoms. 
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the first ionization potential of the Ne atom the DOS's 
show positive values. This implies that the Ceo valence 
electrons transfer resonance energy to the photoelectron 
through intershell coupling. Otherwise the resonance of 
the Ceo is perturbed and the DOS's are reduced by the 
Ne atom. 

Figure 6 shows that when the photon energy is large 
the DOS's are positive; otherwise the DOS's of the Ceo 
are reduced by the Ar atom. This causes negative DOS's 
of the Ar atom encaged in the Cgo ■ 

The solid curves in figures 2, 5 and 6 demonstrate that 
the DOS's of the atoms encaged in the Ceo are positive 
only when the photon energies are large. This indicates 
that strong energy transfer can occur only between the 
Ceo valence electrons and the photoelectron. 

IV. CONCLUSION 

In conclusion, the absolute DOS's for the photoabsorp- 
tion of the Ne, Ar, and Xe atoms encapsulated in the Ceo 
have been calculated in the energy regions both inside 
and outside the Cgo giant resonance. Within the Ceo gi- 
ant resonance the spectra of the Nc, Ar, and Xe atoms 
encaged in the Ceo demonstrate strong oscillation with 
significant increase of the DOS's. This large enhance- 
ment of the DOS's cannot be explained by confinement 
resonance but is due to the energy transfer from the Ceo 
valence electrons to the photoelectron through the inter- 
shell coupling. The energy transfer occurs only when the 
photon energy is large enough to ionize the atom encap- 
sulated in the Cgo- To study the confinement resonance 
it is recommended to avoid the Ceo resonance region by 
selecting the photon energy larger than 40 eV. 
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FIG. 6. (Color online) The symbols have the same meaning 
as in figure 5, except that the curves represent the data for 
the Ar atom. 



of the Ceo giant resonance. The solid curve in Fig. 5 can 
be divided roughly into two parts, a positive part and a 
negative part. When the photon energy is larger than 



This work was supported by the U.S. DOE, Division 
of Chemical Sciences, Geosciences and Biosciences, Of- 
fice of Basic Energy Sciences, Office of Energy Research, 
AFOSR and Army Research Office (Grant W911NF- 
11-1-0194). Calculations used Kraken System (account 
number TG-DMR1 10034) of the National Institute for 
Computational Science, The University of Tennessee. 



[1] Zhifan Chen and A. Z. Msezane, Eur. Phys. J. D 66, 184 
(2012). 

[2] Zhifan Chen and A. Z. Msezane, J. Phys. B submitted 
(2012). 

[3] DM0I3, Accelrys Software Inc., San Diego, CA (2010). 
[4] R. M. Martin, Electronic Structure: Basic Theory and 

Practical Methods, (Cambridge University Press, UK 

2004). 

[5] P. Giannozzi, et al, J. Phys.: Condens. Matter 21, 
395502 (2009). 



[6] B. Walker, A. M. Saitta, R. Gebauer, and S. Baroni, 

Phys. Rev. Lett. 96, 113001 (2006). 
[7] D. Rocca, R. Gebauer, Y. Saad, and S. Baroni, J. Chem. 

Phys. 128, 154105 (2008). 
[8] S. Baroni, et al J. Phys.:Condens. Matter 22, 074204 

(2010). 

[9] M. J. Puska and R. M. Nieminen, Phys. Rev. A 47, 1181 
(1993). 

[10] M. Ya. Amusia, A. S. Baltenkov, L. V. Chernysheva, Z. 
Felfli and A. Z. Msezane, J. Phys. B 38, L169 (2005). 



5 



[11] M. Ya. Amusia, L. V. Chernysheva, and V. K. Dolmatov, 

Phys. Rev. A 84, 063201 (2011). 
[12] M. E. Madjet, et al Phys. Rev. A 81, 013202 (2010). 
[13] Zhifan Chen and A. Z. Msezane, Eur. Phys. J. D 65, 353 

(2011). 

[14] A. L. D. Kilcoyne, et al Phys. Rev. Lett. 105, 213001 
(2010). 

[15] M. Ya. Amusia and A. S. Baltenkov, Phys. Rev. A 73, 
062723 (2006). 

[16] M. E. Madjet, et al Phys. Rev. Lett. 99, 243003 (2007). 
[17] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 78, 1396 (1997). 



[18] A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. 

Joannopoulos, Phys. Rev. B 41, 1227 (1990). 
[19] E. Westin, et al J. Phys. B:At. Mol. Opt.Phys. 29, 5087 

(1996). 

[20] J. Berkowitz, J. Chem. Phys. Ill, 1446 (1999). 
[21] J. Kou, et al Chem. Phys. Lett. 374, 1 (2003). 
[22] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 
[23] W. F. Chan, G. Cooper, X. Guo and C. E. Brion, Phys. 
Rev. A 45, 1420; 46, 149 (1992). 



